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28 A Titan-like organic haze has been hypothesized for Earth’s atmosphere prior to 

29 widespread surface oxygenation ~2.45 billion years ago (Ga). We present a high-resolution 

30 record of quadruple sulfur isotopes, carbon isotopes, and Fe speciation from the -2.65-2.5 

3 1 Ga Ghaap Group, South Africa, which suggest a linkage between organic haze and the 

32 biogeochemical cycling of carbon, sulfur, oxygen, and iron on the Archean Earth. These 

33 sediments provide evidence for oxygen production in microbial mats and localized 

34 oxygenation of surface waters. However, this oxygen production occurred under a reduced 

35 atmosphere which existed in multiple distinct redox states that correlate to changes in 

36 carbon and sulfur isotopes. The data are corroborated by photochemical model results 

37 that suggest bi-stable transitions between organic haze and haze-free atmospheric 

38 conditions in the Archean. These geochemical correlations also extend to other datasets, 

39 indicating that variations in the character of anomalous sulfur fractionation could provide 

40 insight into the role of carbon-bearing species in the reducing Archean atmosphere. 

41 The Campbellrand-Malmani carbonate platform of the -2.65-2.5 billion years old (Ga) 

42 Ghaap Group (Transvaal Supergroup, South Africa 1 ) is one of the oldest carbonate platforms on 

43 Earth. Molybdenum concentrations and isotopes". Re/Mo ratios’ and N isotope systematics" all 

44 suggest biological production of CE in surface waters during deposition of the upper part of the 

45 succession at -2.5 Ga. This suggestion of significant Oz production prior to Earth’s first Great 

46 Oxidation Event ' is supported by nutrient and trace metal systematics in coeval successions from 

47 other areas 6 ' 9 . However, production of O 2 apparently did not lead to a pervasive rise in 

48 atmospheric oxygen levels . 

49 Here, we focus on the lower part of the Ghaap Group, spanning the top of the Boomplaas 

50 Formation through the lower part of the Upper Nauga Formation (Fig. 1 ). Our samples come 
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51 from drillcore GKF01, and represent deposition primarily below wave base in a slope 

52 environment, with cyclic units of microbialites passing upwards into slope carbonates'. 

53 Radiometric age constraints place the base of the core at -2.65 ± 0.08 Ga and the top of the 

54 measured section (within the Upper Nauga) at ~2.5 Ga U! . The succession is exceptionally well- 

55 preserved, having experienced only gentie tectonic warping and sub-greenschist facies 

56 metamorphism ! 1 . 

57 Oxygenation of the Neoarchean oceans 

58 Fe speciation analyses reflect water column redox conditions at the time of deposition, 

59 and show significant variations throughout the section (Fig. 1 ; Table SI). Many samples have 

60 elevated ratios of highly reactive Fe (Fe carbC ma t es + Fe ox , ()es + Femagnetite + Fe^nte) to total Fe 

61 (Fchr/Fct > 0.38), indicating deposition from an anoxic water column 1 ". For those samples 

62 exhibiting high Fchr/Fct, the ratios of pyrite to highly reactive Fe (Fepy/FenjO generally fall 

63 below 0.7-0.8, suggesting that when the bottom waters were anoxic, ferruginous (Fe(II)-rich) 

64 rather than euxinic (sulfide-rich) conditions were dominant '. Other samples representing 

65 shallow water deposition (microbialites of the Monteville Fm) or transported shallow-water 

66 deposits (the Lokammona), fall close to the average oxic Phanerozoic signal 12 , suggesting 

67 formation from an oxic water column" . The low FeHR/Fer ratios and correlative facies 

68 assemblages imply that cyanobacterial 0 2 production in microbial mats resulted in oxygenation 

69 of shallow surface waters as early as -2.65 Ga. 

70 Cyanobacterial oxygen production is further supported by trends in 5 34 S and A 33 S. Sulfur 

7 1 isotope values for the majority of the section fall on or near a linear trend in 5 34 S versus A 33 S that 

72 has been described in Arehean samples (A 33 S = 0.9 x 8 ,4 S)* 4 (Fig. 2A). This trend is interpreted 

73 to reflect mixing between sulfides formed from two primary' atmospheric sources carrying 
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74 opposite sulfur mass independent fractionation (S-MIF) signals - e.g., an oxidized source 

75 (sulfate) with negative A 33 S and a zero valent source (Sg) with positive A 33 S. Samples associated 

76 with microbial mats of the upper Monteville and the Lower Nauga have positive A J, S values that 

77 plot parallel to the primary Archean array, but are enriched in 34 S by ~3 to 796o (Fig. 2A). This 

78 trend could not result from changes in sulfate reduction fractionations alone, as this scenario 

79 would only rotate the array. We suggest this association reflects formation of pyrite from a 

80 residual pool of sulfide generated within the microbial mats’' that was partially oxidized via O 2 

8 1 near the mat surface, enriching the sulfide in 3<i S (see Supplementary Information, SI, for further 

82 discussion). 

83 Redox state of the Neoarchean atmosphere 

84 Covariations in A 33 S, A 36 S, and carbon isotope data reveal an important connection 

85 between the Neoarchean sulfur and carbon cycles. Quadruple sulfur isotope values vary 

86 considerably upcore (Fig. 1 ; Table S2), but all show S-MIF signals reflecting photochemical 

87 production' 0 in an atmosphere devoid of significant 0 2 ' . Samples from the Boomplaas and 

88 Lokammona Formations form a linear array in A 33 S versus A 36 S, with a slope in A 36 S/A j3 S of 

89 approximately -0.96 (Fig. 2B). This slope is close to the reference line defined by previous 

90 measurements of Archean sulfide and sulfate minerals ( — 0.9), and is also inteipreted to reflect 

91 mixing of sulfur from two or more atmospheric sources 16 . In several distinct intervals within the 

92 overlying section, the samples preserve A 33 S and A j6 S values that form arrays rotating in a 

93 clockwise manner, forming steeper slopes in A ,6 S/A 33 S, approaching - 1 .5 (Fig. 2B). The data 

94 preserve a continuum of A 36 S/A j3 S slopes between ~-0.9 and - 1 .5, and extend to highly positive 

95 A 33 S regardless of A 36 S/A 33 S. Biological sulfur cycling can produce scatter in the A 36 S intercept, 

96 since sulfate reduction can affect A 3t 'S much more than A 33 S ! * (Fig. S 1 ). However, a rotation of 
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97 the arrays to different A 36 S/A 33 S slopes instead implies a change in the nature of the primary S- 

98 MIF signature, either due to a change in A 36 S/A 33 S for a single atmospheric signal, or due to 

99 variable amounts of mixing between two or more atmospheric sources with different A 36 S/A 33 S. 

100 Samples w'ith the greatest deviation in A 3h S/A 33 S occur in three distinct zones, at around 1400 m, 

101 1 100 m, and 800 m depth in the core. These zones correspond to intervals exhibiting highly 13 C- 

1 02 depleted organic matter, falling to 5 13 Corg values low'er than -40%o (Fig. 1 ; Fig. S2). Negative 

103 excursions in 8 13 Corg are typically interpreted to reflect higher CFL) fluxes and an increased 

1 04 contribution of methanotrophic biomass to sedimentary organic matter. A significant amount of 

105 methane oxidation (methanotrophy) in the modem oceans occurs anaerobically ly ; however, 

106 FenR/Tex ratios near Phanerozoic oxic values in some of the !3 C-depleted samples suggest that 

107 methanotrophy could also have proceeded via C) 2 (Fig. 1). 

108 A closer examination of previously published Neoarchean sulfur isotope data reveals 

109 similar correlations between changes in A 36 S/A 33 S slope and 8 13 Corg in -2.5 Ga sections in 

1 1 0 South Africa (the Gamohaan Formation) and in Western Australia (the Mount McRae Shale) 

1 1 1 (Fig. S3). Samples forming an array with A 36 S/A 33 S of -1 .5 have also been reported for the -2.73 

112 Ga Tumbiana Formation of the Fortescue Group of Western Australia which preserves a large 

1 1 3 global negative 8 1 J Corg excursion down to -60%o 21 . We suggest that the correlation between 

1 14 negative 8 13 Corg excursions and A 36 S/A 33 S anomalies in two sections and over -200 million 

1 1 5 years implies a global connection between methane and variations in atmospheric S-MIF signals. 

1 16 Modeling the Neoarchean atmosphere and hypotheses for S-MIF production 

1 1 7 Interpreting the magnitude of S-MIF in the geologic record requires knowledge of the 

! 18 production, transformation, and mass partitioning of atmospheric sulfur, as well as subsequent 
1 19 biogeochemical pathways to sedimentation. Recently measured^ absorption spectra for 32 SOj 
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] 20 33 S0 2 , and 34 S0 2 help constrain S-MIF production by S0 2 photolysis and imply that changing 

1 2 i concentrations of UV absorbers'" (such as organic haze 20 "" formed at high atmospheric 

1 22 CH 4 :C0 2 ) might influence the S-MIF signature. S-MIF from S0 2 photolysis could also be 

123 supplemented by other mass-independent kinetic isotope effects " , in reactions such as SO 

124 photolysis, S+S 2 — > S 3 , CS 2 photolysis" 8 "', or the production/destruction of excited-state S0 2 . 

125 Generation of significant S-MIF via any reaction other than SO 2 photolysis remains hypothetical 

1 26 but worthy of future research given that multiple source reactions might help explain the data". 

127 S-MIF, once created, partitions into all atmospheric sulfur species, and requires the simultaneous 

1 28 presence of two (or more) atmospheric exit channels for preservation . Furthermore, these exit 

1 29 channels vary' as a function of atmospheric methane, oxygen, and sulfur concentrations 1430 , and 

1 30 biogeochemical processes could further transform the signal prior to preservation in sedimentary 

131 environments 

1 32 To evaluate these possible alternatives, we ran 1 -D early Earth photochemical models 

133 where we increased CH 4 concentrations until an organic haze formed. The model utilized here 

1 34 advances similar studies" 4 "" by analyzing hazy atmospheres after the evolution of oxygenic 

1 35 photosynthesis. Furthermore, we consider the effect of fractal aggregate scattering by 

136 hydrocarbon aerosols \ and contrast with traditional Mie scattering results. We use these models 

1 37 to explore how optical depth, reaction pathways, and sulfur exit channels vary as a function of 

1 38 CH 4 concentration, and tie these atmospheric properties to candidate mechanisms for the 

139 observed changes in A 36 S/A J, S. Figures 3 and S 8 show two sets of photochemical model 

140 simulations, both varying CH 4 concentrations from 1 to 50,000 ppm, while holding C0 2 

141 concentrations at 1 % and ground-level 0 2 concentrations at 10 ppb. The two simulations differ 

142 only in their hydrocarbon particle scattering physics (see SI for a discussion of boundary' 
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143 conditions, particle formation and scattering). For both fractal and Mie scattering cases, higher 

144 CH4 resulted in increasingly reduced atmospheres, as evidenced by decreasing OH 

145 concentrations (Fig. 3A; S 8 A). Despite large volcanic SO?, and H 2 S fluxes, none of our model 

146 atmospheres accumulated large enough concentrations of SO 2 or OCS to affect S-MIF (see SI for 

147 additional discussion). The dominant contributors to opacity between 180 and 220 nm in our 

148 haze-free models were C0 2 absorption and Rayleigh scattering. Prior to haze formation, most 

1 49 photons longward of 200 nm reached the Earth surface. 

1 50 Even when CH 4 is too low to promote haze formation (CH-uCCF < 0. 1), changes to CH 4 

1 5 1 concentrations altered the redox chemistry and impacted S-MIF exit channels. Enhanced CH 4 

1 52 concentrations alone do not affect S-MIF, as the CH4 absorption spectrum does not overlap the 

153 vibrational bands of SO 2 that give rise to S-MIF. Instead, increases in CH 4 provided the reducing 

1 54 capacity to polymerize elemental sulfur, thus promoting aerosol Sg as an exit channel. 

155 Simultaneously, SO 2 replaced sulfate as the most oxidized exit channel (Fig. 3C; S 8 C), 

1 56 implicating the SCh/Sg couple as carriers of opposite-sign S-MIF. Reaction rates for S 3 

1 57 polymerization, excited state S0 2 formation and SO photolysis all increased with increasing CH 4 

1 58 concentrations (Fig. 3B; S 8 B). If any of these reactions produces significant S-MIF, they could 

1 59 modify the signal produced by SO 2 photolysis. Regardless of mechanism, we suggest this haze 

160 free regime is the origin of the 0.9 Archean ‘reference’ slope in A 36 S/A 33 S, 

161 At CFl 4 :C 02 > 0.1 , a hydrocarbon haze dramatically altered our model atmospheres. The 

162 broadband absorption by this haze caused a drastic reduction in photochemical reactions, and 

163 specifically a slowdown of the S cycle. In the fractal scattering case, SCO photolysis rates slowed 

164 by nearly two orders of magnitude (Fig 3B), as photons shortward of 220 nm were scattered by 

165 the stratospheric haze. By contrast, SO 2 photoexcitation to ‘802(210-327 nm) and 3 S0 2 (337- 
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166 385 nm) remained prevalent. There are two potential sources of S-M1F variation in this regime: a 

1 67 changing ratio of SO 2 photolysis to photoexcitation, and opacity effects arising from fewer total 

1 68 photons. Other pathways involving multiple S atoms, such as S 3 polymerization or SO-SO dimer 

1 69 formation, became unfeasible in the photon-starved CH^COj > 0.1 environment. 

1 70 For 0.1 < CH 4 :C 02 < 0.2, haze-induced shielding led to counter-intuitive decreases in the 

171 CH 4 fluxes needed to support increasing CH 4 concentrations. As atmospheric CH 4 rose the haze 

172 thickened, decreasing both CH 4 photolysis and the formation of radicals that destroy CH 4 . The 

1 73 haze therefore diminished the photochemical sinks for CH 4 such that larger concentrations could 

174 be maintained at lower CB 4 source fluxes. This would have made atmospheres with 0.1 < 

1 75 CFLgCCb <0.2 unstable. If a biosphere produced enough CH 4 to reach CITgCOi > 0. 1 , a positive 

1 76 feedback would have ensued, increasing CH 4 concentrations until CH 4 :CC >2 = ~ 0.2, at which 

1 77 point further increases in CH 4 concentrations would have required biologically unfeasible CFI 4 

1 78 fluxes (see details in SI). 

1 79 In both simulations, atmospheres with CH 4 .CO 2 > 0.2 are unsustainable. If organic 

1 80 particles scatter light according to Mie theory, a thick haze could be prevented by a Gaian 

181 climate/CHU feedback*^ on global glaciations . If organic particles behave like fractal scatterers, 

182 a thick haze would have formed a UV shield effective enough to shut down photochemical haze 

1 83 production (see Fig. 3B and SI). In either case, we can place an upper limit of ~0.2 on the 

184 atmospheric CH^rCC^ ratio. 

1 85 In summary, there are two accessible, stable regimes in our photochemical simulations: 

186 CH 4 .CO 2 <0.1, and CHL^COj ~ 0.2. In the first of these regimes, SO 2 photolysis was a major 

1 87 source of S-MIF, although other potential S-M1F reactions proceeded rapidly enough to be 

1 88 quantitatively relevant. The exit channels varied as a function of atmospheric redox state, and the 
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1 89 magnitude of S-MIF in a given exit species would be expected to vary inversely with the relative 

190 mass fraction. In the second regime, the atmospheric S cycle operated slowly due to haze- 

1 9 1 induced opacity that significantly decreased photolysis rates in the lower atmosphere. S0 2 

1 92 photolysis remained a likely source of S-MIF, and the lower photon availability could have led to 

193 significantS02 isotopologue absorption effects. S-MIF from symmetry-dependent de-excitation 

1 94 of S0 2 also had the potential to be an important source of S-MIF in this regime. S0 2 became the 

195 dominant exit channel, so its S-MIF signature should trend to 0 by mass balance considerations, 

196 with correspondingly higher S-MIF magnitudes in the minor (by mass) exit channels. Transition 

1 97 between the first and second regimes could therefore induce changes in S-MIF source 

1 98 magnitudes, coupled with changes in the reservoirs by which S-MIF exits the atmosphere. If low 

1 99 b' Corg reflects periods of enhanced methane flux, a bi-stable atmospheric transition from haze- 

200 free (CH 4 .CO 2 < 0. 1 ) to thin haze (CH4:C0 2 = 0.2) could explain the A 36 S/A 33 S slope changes in 

201 ancient sediments. 

202 Consideration of these models along with the high-resolution geochemical dataset 

203 presented here points to a Neoarchean depositional environment that includes: (1) microbial mat 

204 ecosystems similar to modem mats, with intensive internal recycling of sulfur; (2) cyanobacterial 

205 oxygen production and shallow water oxygenation as early as -2.65 Ga; and (3) a link between 

206 methanotrophy and S-MIF producing atmospheric chemistry, indicating a reducing atmosphere 

207 with multiple episodes of organic haze formation during periods of enhanced biological CFU 

208 production. These records provide crucial evidence for the biogeochemical production and fate 

209 of oxidants and reductants in the oceans and atmosphere during the critical period in Earth 

210 history' immediately pre-dating the first major rise in atmospheric oxygen. 

211 
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212 Methods 

213 Pyrite sulfur isotope compositions were determined on Ag 2 S precipitates from distillation 

2 1 4 with chromous chloride. Quadruple S isotope measurements were performed at the University of 

215 Maryland, following the techniques described in Domagal-Goldman et air 4 , and reported as S 34 S 

216 = (( 34 S/ 32 S) sanip ,e/( 34 S/ 32 S)v-cDT - 1 ), A 33 s = ( 33 S/ 32 S) sampi e/( 33 S/ 32 S)v.cDT-[( 34 S/ 32 S) sampie /( 34 S/ 32 S)v. 

217 corf 515 , A 36 s = ( 36 S/ 32 S) sample /( 36 S/ 32 S)v-cDT-[( 34 S/ 32 S) sampl ,/( 34 S/ 32 S)v-cD T ] 1 9 , with values given in 

2 1 8 permil (96o). Analytical uncertainties on S isotope measurements, estimated from long-term 

219 reproducibility of AgjS fluormations, are 0.14, 0.008, and 0.20 (la) for 6 34 S, A 33 S, and A 36 S, 

220 respectively. Iron speciation was determined at Newcastle University via the sequential 

221 extraction techniques of Poulton et al. l2 u . Replicate extractions give an RSD of <5% for all 

222 extraction steps. Total organic carbon (TOC) was measured on a Leco analyzer after treatment 

223 with dilute 20% HC1 to remove carbonate phases. Organic C isotopes were measured on acid- 

224 treated samples by EA-IRMS, reported as 8 I3 C = (( r, C/ J 2 C) sam ,, : , ; /( 13 C/ 12 C) V pdb - 1) (%o). Replicate 

225 analyses (n = 6) of an internal standard (5 ,3 C org = -26.43%o) gave an average value of -26.43 ± 

226 0.0696o ( 1 a) during the sample run. 

227 The 1-D photochemical model derives from Zhanle et al. "\ but contains the organic 

228 carbon species and chemistry' from Domagal-Goldman et al 2<t , as well as their 5 two-stream 

229 radiative transfer scheme. The mode! simultaneously solves photochemical production and loss 

230 for 74 species undergoing 392 reactions, including transport by eddy and molecular diffusion 

231 over a 100 km grid with 0.5 km grid spacing. We adopt the temperature, eddy diffusion, and 

232 water vapor profiles from Domagal-Goldman et al , 24 , leaving all other parameterizations 

233 (rainout, lightning, diffusion-limited hydrogen escape, etc.) unchanged. The modem solar flux is 
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234 used, although sensitivity tests were run approximating the spectral character of the sun at 2.5 

235 Ga. Further details, including boundary' conditions, are given in the Supplementary' Information. 

236 

237 Correspondence and requests for materials should be addressed to Aubrey L. Zerkle, 

238 aubrey.zerkle@ncl.ac.uk. 
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Figure Legends 

Figure 1. Lithologic and geochemical data for the bottom part of core GKF0 1 , through -2.65 to 
2.5 Ga sediments of the Ghaap Group (lithologies from Schroder et al. 1 ). Data include pyrite 
sulfur isotope values (§ 34 S, A 33 S, A 36 S, and A 36 S/A j3 Sdev, calculated as the absolute value of the 
deviation from a A 36 S/A 33 S slope defined by the Boompiaas and Lokammona Formations), 
carbon isotopes of organic carbon (8 13 Corg), and Fe speciation data (Tables S2 and S3). Symbols 
are defined in Figure 2. Grey circles are from Ono et airy grey squares are from Kendall et al. 4 , 
and grey diamonds are from Fischer et al. M) . 

Figure 2. Cross-plots for quadruple sulfur isotopes from the GKF01 section. A. S 34 S versus A 33 S 
values. The dotted line is the Archean reference line (A j3 S = 0.89 * S 34 S) 14 . B. A 33 S versus A % S 
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355 values. The solid line represents a reference array with A 36 S/A 33 S slope of -0.9, and the dotted 

356 line represents an array with a slope of -1 .5. 

357 Figure 3. Photochemical model results for fractal particles. Models were run with COj fixed at 

358 1% and ground-level CH 4 mixing ratios from IxlO ' 6 to 0.05, displayed as CH^rCO? ratios. 

359 Shaded regions indicate unconverged solutions. A. Mixing ratios (solid lines, left axis) and 

360 fluxes (dashed lines, right axis). B. Photochemical reaction rates relative to the rate at CFLpCC^ = 

361 I O ' 4 (non-red lines, left axis). Effective optical depth at 200 nm (solid) and 550 nm (dashed) (red 

362 lines, right axis). C. Sulfur exit channels shown as fraction of total sulfur leaving the atmosphere. 

363 The dotted line is the sum of exit channels shown. 
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Figure 3. Zerkle et al. 




